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Abstract 

Our understanding of magnetic reconnection in resistive magnetohydrodynamics has gone 
througli a fundamental change in recent years. The conventional wisdom is that magnetic re- 
connection mediated by resistivity is slow in laminar high Lundquist (S) plasmas, constrained by 
the scaling of the reconnection rate predicted by Sweet-Parker theory. However, recent studies have 
shown that when S exceeds a critical value ~ 10^, the Sweet-Parker current sheet is unstable to 
a super- Alfvenic plasmoid instability, with a linear growth rate that scales as S^^^. In the fully 
developed statistical steady state of two-dimensional resistive magnetohydrodynamic simulations, 
the normalized average reconnection rate is approximately 0.01, nearly independent of S, and the 
distribution function f{ijj) of plasmoid magnetic flux ip follows a power law f{tp) ~ ip~^ . When Hall 
effects are included, the plasmoid instability may trigger onset of Hall reconnection even when the 
conventional criterion for onset is not satisfied. The rich variety of possible reconnection dynamics 
is organized in the framework of a phase diagram. 
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I. INTRODUCTION 



Magnetic reconnection is generally believed to be the underlying mechanism that powers 
explosive events such as flares, substorms, and sawtooth crashes in fusion plasmas. 
Traditionally, magnetic reconnection is grossly classifled into two categories, i.e. collisional 
and collisionless reconnection, depending on whether ions and electrons remain coupled or 
not in the diffusion region. The conventional wisdom is that laminar collisional reconnection 
is described by the classical Sweet-Parker theory, which assumes an elongated current 
sheet characterized by a length L of the order of the system size. The governing dimensionless 
parameter is the Lundquist number S = VaL/t], where Va is the upstream Alfven speed, L 
is the reconnection layer length, and rj is the resistivity. According to Sweet-Parker theory, 
the reconnection rate ~ BVa/ VS, where B is the upstream magnetic fleld. In many systems 
of interest, the Lundquist number S is very high (e.g. S ~ 10^^ — 10^^ in solar corona), 
and the corresponding Sweet-Parker reconnection rate is too slow to account for energy 
release events. For this reason, research on fast reconnection in the past two decades has 
mostly focused on collisionless reconnection, which can yield reconnection rates as fast as 

The Sweet-Parker theory assumes that the elongated current sheet is stable. However, 
it has long been known that this assumption may not hold, because the current sheet can 
become unstable to tearing instability and spontaneously form plasmoids at high-S*. jo - 12| 
Biskamp estimated the critical Lundquist number Sc to be approximately lO'^.jol Obser- 



vational evidences o 
solar atmosphere. Il3 



plasmoids have been reported in the Earth's magnetosphere and the 



- |23] Although numerical and observational evidences of plasmoids were 
abundant, precise scalings of of the linear instability (hereafter the plasmoid instability) 
were not known until recently. In a seminal paper, Loureiro et al. predicted that the linear 
growth rate 7 scales as 7 ~ S^^^L/Va, and the number of plasmoids scales as S''^/^.|24| At 
flrst sight, these results are rather counterintuitive, because linear growth rates of most re- 
sistive instabilities scale with 5* to some negative fractional power indices instead of positive 
ones. The crucial point is to realize that the equilibrium, i.e. the Sweet-Parker current 
sheet, also scales with 5*. The results of Loureiro et al. can be readily derived once the 
Sweet-Parker scaling of the current sheet width 6sp ^ L/\/S is incorporated in the classical 



tearing mode dispersion relation. 



25| 



While the work of Loureiro et al. drew attention to the surprising scahng properties of 



the hnear plasmoid instabihty, 



24| it is the nonhnear behavior of the instabihty that has 



proved to be transformational for traditional reconnection theory. The onset of the linear 
plasmoid instability would not be nearly as interesting were it not for the fact that it leads to 
a new nonlinear regime, 
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29| entirely unanticipated by Sweet-Parker reconnection theory. 



Furthermore, the plasmoid instability, which appears to be ubiquitous in high-S* plasmas, 
often leads inevitably to kinetic regimes in which the onset of fast reconnection occurs earlier 
than was previously thought possible. 
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32| The picture emerging from these recent studies 



is that large-scale, high-Lundquist-number magnetic reconnection can exhibit complex and 
rich dynamics that we are only beginning to understand. 



The main objective of this paper is to give our perspective of the recent advances in 
this subject, and address some issues that have arisen during the course of this work but 
have not been addressed in previous publications. Most of our discussion is limited to 
two-dimensional (2D) systems. The plasmoid instability in three-dimensional (3D) systems, 
which is currently an active area of research, will be briefly discussed. This paper is organized 
as follows. Section [Til reviews the linear theory, with emphasis on the convective nature of 
the instability that is one of its key features. Section IIIII discusses plasmoid dynamics in 
the fully nonlinear regime. Results on reconnection rate and scaling laws are reviewed, and 
a heuristic argument is given. Section IIVI discusses statistical descriptions of the plasmoid 
distribution, which has been a topic of considerable interest and debate in recent years. We 
approach this problem with a combination of analytical models of plasmoid kinetics and 
direct numerical simulations (DNS), where the plasmoid distribution function is found to 
obey a power law of index —1 for smaller plasmoids, followed by an exponential falloff for 
large plasmoids. The condition for transition from the power-law regime to the exponential 
tail is discussed in great detail, which is crucial for interpreting the numerical results. Section 
IVl addresses the role of the plasmoid instability on the onset of collisionless reconnection. A 

and 



revised phase diagram is presented that explicitly includes the physics of bistability 
the intermediate regime reported in Ref. |32|. Although these effects were discussed before, 
they were omitted in our previous rendition of the phase diagram. Finally, open questions 
and future challenges are summarized in Sec. IVI[ 



II. LINEAR THEORY OF THE PLASMOID INSTABILITY 



Consider a Harris sheet of width a with the equihbrium magnetic field B = tanh(2;/a)x. 
According to the classical linear tearing instability theory, the tearing mode growth rate for 



the "constant--?/^" and "non- constant-'?/'" regimes are as follows, respectively: 



7 ~ — X S , (1) 



a 



where Sa = clVa/tj is the Lundquist number based on the current sheet width a. The 
transition from "constant-?/'" to "non-constant-?/^" modes occurs at kmaxO- ~ where 
the linear growth rate peaks and scales as 7max ~ S~^^'^{VA/a)- Note that the linear growth 
rate is proportional to Sa raised to a negative fractional exponent in both regimes. To apply 
the tearing mode theory to a Sweet-Parker current sheet, it is important to note that the 
current sheet length L is dictated by the global geometry, and remains approximately the 
same when rj varies. On the other hand, the current sheet width 6sp follows the scaling 
6sp ~ L/y/S, where the Lundquist number 5* = LVa/^ is now defined with the current 
sheet length. The current sheet width a in should now be replaced by 5sPi and Sa is related 
to 5* via the relation Sa = ^spVa/ti = S^^"^. After some algebra, the dispersion relation ([!]) 
can be rewritten as 

7 ~ — X < , (2) 

L \ ^y^, K « S'/' 

where k = kL and e = ^sp/L ~ S^^^"^. The peak growth rate occurs at n^ax ~ S^^^, 
with 7max ~ S^/'^{Va/ L). The number of plasmoids generated in the linear regime can be 
estimated as ~ L/\max ~ f^max ~ S*"^/^, where Xmax is the wavelength of the fastest 



„ .node. These se^.„g .elat.ons have been confl.n,ed by nun,e.ic^ studies.bdM^ 

The plasmoid instability is referred to as a super- Alfvenic instability, [25( because the linear 
growth rate far exceeds the inverse of the Alfven time scale Ta = L/Va along the current 
sheet in the high-S* limit. It should be noted, however, that the growth rate is always 
slower than the inverse of the Alfven time scale 5sp/Va transverse to the current sheet, as 
is required for the tearing mode analysis. 



The above analysis (also the analysis in Ref. 2J|) ignores the effects of flow, which is 



self-consistently generated in a Sweet-Parker current sheet, on the instability. However, it is 



important to appreciate the convective nature of the instabihty. Because the outflow speed 
is of the order of Va, disturbance in the current sheet will be convected out at the Alfven 
time scale r^. Based on the scaling 'jmax ~ S^^^C^a/ L), the fastest growing mode will be 
amplifled by a factor of ~ exp(S'^/^) before being convected out of the layer. If the initial 
disturbance is sufficiently small, it may not grow to a perceptible size, and the current sheet 
will appear to be stable. For this reason, we do not expect a clear-cut critical value of 5* 
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El 



for the instability. Simulations typically give a critical value Sc ~ 10 
although clean simulations that remain stable up to S* ~ 10^ have been reported. 138| There 

n 

are also indications that the critical value Sc may depend on plasma /3.[39j Because the linear 
amplification factor increases rapidly as S becomes large (e.g., exp{S^/^) ~ 10^ 10^^, 10^^ 
at 5* = 10^,10^,10^, respectively), the Sweet-Parker current sheet becomes very fragile at 
high S, and plasmoid formation, which is the consequence of a true physical instability, is 
unavoidable. For instance, while noise levels significantly lower than 10"^ of the background 
field may be sufficient to obtain a stable Sweet-Parker current sheet at 5* = 10^, it requires 
a noise level significantly lower than 10"^^ to achieve a similar level of stability at S* = 10^. 
An even more stringent upper bound on the noise level is obtained at higher values of S. 
We should point out that the linear amplification factor is introduced here to estimate the 
requirement on noise in order to keep the Sweet-Parker current sheet stable. It may not be 
representative of the real amplification, because the initial exponential growth is expected 
to slow down when nonlinear effects become important, which typically occurs once the the 
plasmoid size becomes comparable to the current sheet width. (40| 



III. EFFECTS OF PLASMOID INSTABILITY IN COLLISIONAL MAGNETIC RE- 
CONNECTION 



After onset of the plasmoid instability, the reconnection layer changes to a chain of 
plasmoids connected by secondary current sheets that, in turn, may become unstable again. 
This process of cascading to smaller scales is reminiscent of fractals. {ll| If the large-scale 
configuration evolves slowly, eventually the reconnection layer will tend to a statistical steady 
state characterized by a hierarchical structure of plasmoids (see Fig. [T] for snapshots of 
the cascade). Two-dimensional numerical simulations show that reconnection rate becomes 



nearly independent of S in this regime, with a value ~ O.OlV^i?. 
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37| The number of 
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Figure 1. (Color online) Out-of-plane current density at different times shows fractal- like cascade 
to smaller scales via the plasmoid instability. First, the Sweet-Parker current sheet breaks up to 
form a chain of plamoids connected by secondary current sheets (top panel). Secondary current 
sheets are Sweet-Parker like, become unstable again, and generate the next batch of plasmoids 
(middle panel). This cascade leads to a hierarchy of plasmoids of various sizes (bottom panel). 
The. .ap..o.. p..ese,« a e.all po.io^ „, »ho,e ...latio,. fro. a S ^ 10' ..nu.a.io. 
The reader is referred to Fig. 4 of Ref. [28| for an illustration of the whole system. 



plasmoids rip, the widths 5 and lengths / of secondary current sheets follow scaling relations 
Up (X S, 6 (x 1/S, and / oc 1/S'.|28| These scaling relations may be understood by assuming 
that all secondary current sheets are close to marginal stability. The rationale behind the 
assumption is as follows. Firstly, we note that cascade to smaller scales will stop if the local 
Lundquist number IVa/t] of a secondary current sheet is smaller than Sc- Secondly, secondary 
current sheets typically get stretched and become longer over time due to gradient in the 
outflow, which on average increases from zero near the center to ~ Va_ near the ends of the 
reconnection layer. Thirdly, when the current sheet length / becomes sufficiently long such 
that IVa/ti > Sc, the current sheet becomes unstable, new plasmoids are generated, and the 
fragmented current sheets become short again. Consequently, we expect the local Lundquist 
number IVa/t] to stay close to Sc, i.e. I ~ tiSc/Va ~ LSc/S. The corresponding current sheet 

6 



width S ~ 1/ ^fSl ~ LSy^/S, and the number of plasmoids is estimated asrip L/l S/Sc- 
Finally, reconnection rate can be estimated by r]J r]B/5 BVA/y/Sl ~ IO'Va^, which 
is independent of S. These scaling relations are consistent with the simulation results. 



IV. STATISTICAL DISTRIBUTION OF PLASMOIDS 



Statistical descriptions of plasmoids have drawn considerable interest in recent years 
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42] partly due to the possible link between plasmoids and energetic particles. 
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43| The 



scales, 



ractal-like fragmentation of current sheets suggests self-similarity across different 



ll| which often gives rise to power laws.j4^ A recent heuristic argument by Uzdensky 



et al. suggests that if we consider the statistical distribution of the plasmoids in terms of 
their magnetic fluxes the distribution function f{ip) follows a power law f{'4)) ~ ip''^. 29| 
Th,s result canbe formally derived by adopting a mode, of plasmo.d U.neties (similar tothat 



given in Ref. 



4l| ) and obtaining steady-state solutions of the plasmoid distribution. 

m 



The governing kinetic equation for the time evolution of f{'4)) is written as 

7^ + "77T = C^W , (3) 

ot dip Ta Ta 

where N{%Ij) = f{ip')dip' is the cumulative distribution function, i.e. the number of 
plasmoids with fluxes larger than ip. Several idealized assumptions have been made in 
writing Eq. ([3]). Firstly, the flux of a plasmoid grows due to reconnection in adjacent 
secondary current sheets. Following the assumption that all secondary current sheets are 
close to marginal stability, the flux of a plasmoid grows approximately at a constant rate a ~ 
-BVa/a/S^- This gives the plasmoid growth term ad f /dtp on the left hand side. Secondly, 
new plasmoids are created when a secondary current sheet becomes longer than the critical 
length for marginal stability. We assume that when new plasmoids are created, they contain 
zero flux; this is represented by the source term C'^('^)) where 6{ip) is the Dirac 5-function, 
and C is the magnitude of the source. This source term sets the boundary condition for 
f{ip) at ip = 0. Thirdly, plasmoids disappear due to coalescence with larger plasmoids, 
which is assumed to be instantaneous. Assuming the characteristic relative velocity between 
plasmoids is of the order of Va, the time scale of a plasmoid with flux to encounter a 
larger plasmoid is estimated as ~ L/N{iI))Va ~ Ta/N{iIj). This gives the coalescence loss 
term —JN/ta- Note that when two plasmoids coalesce, the flux of the merged plasmoid 



is equal to the larger of the two original fluxes. |4l| Therefore, coalescence does not affect 
the value of / at the larger of the two fluxes. Lastly, plasmoids are advected out from the 
reconnection layer with speeds ~ on a characteristic time scale ta- This is represented 
by the advection loss term — //r^. 



Exact steady-state solutions of Eq. ([3]) can be found analytically. |45| However, for the 
discussion here it is instructive to consider approximate solutions in different regimes. At 
large ip when ^ 1, the steady-state equation reduces to ad f /Oil) ~ —J/ta- In this regime 
/ ~ exp(— ?/'/arA). On the other hand, when ^ 1, the advection loss term is negligible, 
and we have ad f /dip ~ —JN/ta- In this regime, A^ ~ 2aTA'ip~^ and / = —ON /dip ~ 
2aTAip~'^ is the solution. As such, the steady state solution admits both an exponential tail 
and a / ~ ip""^ power-law regime. The dominant loss mechanism in the former regime is 
advection, while it is coalescence in the latter. In other words, the plasmoids in the power- 
law regime must be deep in the hierarchy, whereas large plasmoids follow a distribution that 
falls off exponentially. Transition from the power-law regime to the exponential tail occurs 
when A^ ^ 1, i.e. at ■?/; ^ 2aTA- The distribution function f{ip) also deviates from the 
ip'"^ power law in the small ip limit, otherwise the cumulative distribution function N{ip) 
will diverge as 0. Because N{ip) Up S/ Sc &s ip ^ 0, the transition occurs 

when 2aTAip~'^ — np, i.e. when ip ^ 2aTA/np. Therefore, the power law holds in the range 
2aTA/np <^ip <^ 2aTAi which becomes more extended for higher 5*. 

This prediction of f{ip) ~ ip~'^ power-law distribution can be tested with numerical 
simulations. Figure [2] shows the cumulative distribution function N{ip) (panel (a)) and the 
distribution function f{ip) (panel (b)) from a. S = 10^ simulation reported in Ref. 45|. The 
dataset was constituted of 30507 plasmoids collected from an ensemble of 521 snapshots 
during the quasi-steady phase with a cadence of 100 snapshots per ta; an ensemble average 
was carried out for better statistics. The distribution function f{ip) exhibits an extended 
power-law regime; however, the power law is close to f{ip) ~ ip~^ instead of f{ip) ~ ip~'^. 
The vertical dotted line denotes where N{ip) crosses A^ = 1, indicating the switch of the 
dominant loss mechanism from coalescence to advection. From the above analysis, this 
switch of the dominant loss mechanism is responsible for the transition from a power-law 
distribution to an exponential falloff. And indeed, where N{ip) crosses A^ = 1 approximately 
coincides with where the distribution function f{ip) starts to deviate from f{ip) ~ ip~^ to a 
more rapid falloff. Incidentally, this rapidly falling tail was where Loureiro et al. attempted 

8 



10^ 



10' 



_ 



10-1 
10-2 

10-3 
108 

10^ 

10^ 

g 10^ 

10* 

103 

102 

lOi 
10 



*** 



***■ 



(b) 



10- 



10-^ 10-5 



10- 



10- 



10- 



Figure 2. (a) Cumulative distribution function N(ip) and (b) distribution function /(V') of plasmoids 
from a S = 10^ simulation. The vertical dotted line denotes where A^(^) = 1. 



to fit with the f{ip) ~ power law. 37| At smaller ip, their reported distribution also 
appears to be consistent with f{ip) ~ Because our simulation only lasted a few ta, the 
statistics in the large-?/; regime is sufficiently uncertain that it is difficult to make a clear 
distinction between a and an exponential falloff. Observationally, the distribution of the 
flux transfer events (FTEs) in the magnetopause from Cluster data, collected over a period 



42| 



of two years, appears to be consistent with an exponential tail. 

An alternative way to make a distinction between tp''^ and tp'^ distributions is to examine 
the distribution of the leading digit d of the flux. (For example, if ip = 2.35 x 10"^, then 
d = 2.) Although the leading digit surely depends on the unit we use, the distribution of 
the leading digit will remain the same if the underlying distribution is a power law. For the 
f{ip) ~ ip"^ distribution, the probability P{d) of the leading digit d follows Benford's law 
P{d) = logio(l + l/ci).|46| On the other hand, the probability will be P{d) = 10/9d{d + l) if 
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Figure 3. (Color online) The distribution of the leading digit of flux shows good agreement with 
Benford's law (for f{ip) ~ ^~^), but deviates significantly from the prediction based on /{i/j) ~ 



the distribution is /{ip) ~ ip'"^. Figure [3] shows the distribution of the leading digit, which 
is in good agreement with Benford's law, but deviates significantly from the prediction 
corresponding to the distribution f{ip) ~ Because all plasmoids in the dataset, not 

just those in the power-law regime, are used in this analysis, the good agreement with 
Benford's law reflects the fact that the majority of the plasmoids are in the f{ip) ~ ip"'^ 
power-law regime. 

As we have discussed, the dominant balance in Eq. ([3]) leading to the f{ip) ~ solution 
is between the plasmoid growth term and the coalescence loss term, i.e. ad f /dip ~ —JN/ta- 
A key assumption underlying the loss term —JN/ta is that the relative speeds of a plasmoid 
with respect to neighboring plasmoids larger than itself are of the order of Va and are 
uncorrelated to the flux of the plasmoid. However, that was found not to be the case 
when we examined the numerical data. Instead, the relative speeds between large plasmoids 



tend to be lower. 45| Figure H] shows the stack plot of plasmoid positions along the outflow 
direction x over a period of 2tai where each dot represents a plasmoid, color coded according 
to its flux in logarithmic scale. The stack plot has a tree-like structure, with branches that 
represent the trajectories of larger plasmoids shown on the red side of color scale. Those dots 
(or "leaves") on the blue side of the color scale are mostly smaller plasmoids that do not show 
perceivable trajectories. We can clearly see that many branches are nearly parallel to their 
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Figure 4. (Color online) Stack plot of plasmoid positions along the outflow direction x during 
the period t = 4 to t = 6. Each dot represents a plasmoid, color coded according to its flux in 
logarithmic scale. 



neighboring branches, indicating low relative speeds between large neighboring plasmoids. 

This interesting phenomenon may be understood as follows. Roughly speaking, the flux 
of a plasmoid is proportional to its age because all plasmoids approximately grow at the 
same rate a. Consequently, a plasmoid can become large only if it has not encountered 
plasmoids larger than itself for an extended period of time. Those plasmoids moving rapidly 
relative to their neighbors will encounter larger plasmoids and disappear easily, whereas 
those with small relative speeds are more likely to survive for a long time and become large. 
To incorporate this important effect, we have to consider a distribution function not only 
in flux, but also in velocity as well. Let F{iIj,v) be the new distribution function, where v 
is interpreted as the plasmoid velocity relative to the mean flow. We propose the following 
kinetic model for F{iIj,v): 

d,F + a— = CmKv) , 

di) Ta Ta 

where the function H is defined as 



(4) 



dip 



oo I /I 

dv'\^^^:^F{i/j',v'), 



(5) 



and h{v) is an arbitrary distribution function in velocity space when new plasmoids are 
generated. The distribution function f{il)) can be obtained by integrating F['il), v) over 
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Figure 5. (Color online) Top: On the right hand side of the panel is an example of two partially 
merged plasmoids, both small in size. Bottom: When the merging is completed, they become 
a plasmoid much larger in extent. Also there are three new plasmoids generated in the middle. 
Contour lines represent magnetic field lines, and plasmoids are color coded with fluxes in logarithmic 
scale. Note that the two panels cover different ranges in the x direction, because the whole structure 
is moving. 



velocity space. Eq. (jlj) differs from Eq. ([3]) in the coalescence loss term, where the relative 
speed |t; — between two plasmoids is taken into account in the integral operator of Eq. 
([5]). If we replace |w — in Eq. ([5]) by Va, then Eq. (jlj) reduces to Eq. (j3j) after integrating 
over velocity space. By numerically solving for steady-state solutions of Eq. (j4j), we find 
that the distribution in the power-law regime is close to f{ip) ~ ^A"^, consistent with DNS 



results (see Figure 4 of Ref. (45|). This conclusion does not appear to be sensitive to the 
specific form of h{v), as long as h{v) covers a broad range of v (typically of the order of Va)- 

We caution the reader that although the heuristic argument in Sec. IIIII and the kinetic 
model in this Section appear to account for the numerically observed scaling laws and plas- 
moid distributions, they are far from complete and satisfactory descriptions of the complex 
dynamics of a plasmoid-dominated reconnection layer. In particularly, a key building block 
of the present theory is the assumption that plasmoids are connected by marginally stable 
current sheets. That is clearly an oversimplification and does not appear to be always (and 



quite often not) the case upon close examination of simulation data; 28| therefore, further 
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exploration is warranted. 



47| Some other potentially important effects are also omitted. For 



instance, plasmoids are treated as point "particles" and coalescence between them is assumed 
to occur instantaneously, whereas in reality larger plasmoids take longer to merge, and there 
can be bouncing (or sloshing) between them. 48|, |49| Furthermore, the velocity v relative to 
the mean flow is assumed to remain constant throughout the lifetime of a plasmoid, whereas 
in reality some variation is expected due to the complex dynamics between plasmoids. 

The fact that plasmoids can be in a partially merged state can have significant effects 
on the statistical distribution, 37| and this issue is further complicated by that there are 



different ways of identifying the extent of a plasmoid among researchers. In our diagnostics, 
plasmoids are extrema (0-points) of the flux function, which we solve as a primary variable 
in the simulation code. The extent of a plasmoid is determined by expanding the level 
set of the flux function from the extremum until it reaches a saddle point (X-point). Note 
that our convention differs from that adopted by Fermo et al. (see Figure 3 of Ref. 42|). 
Our method treats all plasmoids on equal footing, whereas that of Fermo et al. makes a 
distinction between the dominant and the lesser plasmoids for partially merged plasmoids. 
Although our method is mathematically unambiguous, it has a consequence that when two 
plasmoids are in the process of merging, both shrink in size until the merging is completed 
and a large plasmoid appears suddenly (Figure |5]), whereas in the convention of Fermo et 
al. the lesser plasmoid shrinks in size and the dominant one keeps growing. This is a subtle 
point that merits further consideration if we apply the plasmoid distribution functions to 
the problem of particle energization. 



V. ROLES OF PLASMOID INSTABILITY IN ONSET OF COLLISIONLESS RE- 
CONNECTION 

Thus far, our discussion assumes that resistive MHD remains valid down to the smallest 
scales. This assumption is clearly questionable when current sheet widths reach kinetic 
scales, when two-fluid (Hall) and kinetic effects become important. Conventional wisdom 
had it that the onset of fast reconnection occurs when the Sweet-Parker width 5sp is smaller 
than the ion inertial length di (which should be replaced by the ion Larmor radius at 



the sound speed, p^, if there is a strong guide field). 
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54| Because secondary current 



sheets can be much thinner than the Sweet-Parker width, the implication is that coUisionless 
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Figure 6. (Color online) (a) Phase diagram of magnetic reconnection that shows the five phases, 
(b) Expanded view about S = 10^ and A = 10^ shows two history-dependent regimes. 



'or onset is not met. This 
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reconnection may set in even when the conventional criterion 
has been confirmed by several recent studies using fully kinetic |30[ and Hall MHD 
models. Along with these studies, a practice has emerged of using phase diagrams to describe 
various possible "phases" of reconnection in the parameter space of S and A = L/di. 32|, |55- 
57| Figure |6] shows our current rendition of the phase diagram, which is divided into seven 
regimes (including two history-dependent regimes), and five "phases", that will be detailed 
as follows. 

The collisionless Hall phase is realized when the conventional criterion for onset is satis- 
fied, i.e. 6sp < di. In the parameter space of S and A, that translates io S > h?. This is the 
regime where reconnection dominated by collisionless effects will surely occur. In the region 
S* < A^, coUisional Sweet- Parker reconnection can be realized if the current sheet is stable, 
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i.e. when 5* < Sc- However, that is not always the case. An important phenomenon discov- 
ered by Cassak et al. is that reconnection is not uniquely determined by the parameters S 
and A alone; in a so called "bistable" region of the parameter space, both Sweet-Parker and 
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Hall reconnection are possible, depending on the history in parameter space. 
A useful parameter to discuss this hysteresis phenomenon is the Lundquist number based 
on di, defined as Sd- = VAdi/f] = S/A, which characterizes the effect of resistivity on the di 
scale. Disregarding the plasmoid instability, the condition for bistability may be expressed 
as A > > S^., or equivalently 

ASl <S<A\ (6) 



where is a critical value of 5"^.. The value of S^. is found to be ~ 0(10^); 32l . |33 | 



however, numerical evidence indicates that S^. may increase with increasing A, 32| although 
the precise scaling is not known. The parameter space defined by Eq. ([6]) is represented by 
the shaded region enclosed by dashed lines in Fig. El where we have assumed an arbitrary 
scaling just to show qualitatively the dependence of S^- on A. The regime in the overlapping 
region defined by Eq. ([6]) and S < Sc is the true "bistable" regime, where either collisionless 
Hall or coUisional Sweet-Parker reconnection can be realized, depending on the history (see 
panel (b) of Fig. [6] for a zoom-in view). 

When S > Sc, the Sweet-Parker current sheet is unstable. However, reconnection may 
remain coUisional provided that secondary current sheet widths stay above the di scale. 
Using the scaling estimate 5 ~ LS^^/S, the criterion is 

Sc<S < ASy\ (7) 

which defines the coUisional plasmoid-dominated phase. Note that the there is an overlap- 
ping region defined by Eq. and Eq. ([7]) together (again see panel (b) of Fig. [H] for a 
zoom-in view). That suggests there may exist a history-dependent regime where both col- 
lisionless Hall and coUisional plasmoid-dominated reconnection can be realized. Although 
this regime has not been demonstrated in simulation, it remains an interesting possibility. 

The region defined by AS^^ < S < A^ is where secondary current sheets can reach the di 
scale. Where this region overlaps the "bistable" region defined by Eq. is the "plasmoid- 
induced Hall" phase. Although this regime is supposed to be bistable according to Eq. dH]), 
Sweet-Parker reconnection ceases to exist due to the plasmoid instability, and collisionless 
Hall reconnection is the only possibility. The only difference between the "plasmoid-induced 

15 



Hall" regime and the "collisionless Hall" regime is the process of reaching collisionless Hall 
reconnection. In the former, onset of Hall reconnection is proceeded by cascading to smaller 
scales induced by plasmoids, whereas in the latter the primary current sheet will reach the 
kinetic scales and immediately trigger the onset of Hall reconnection. 



The region defined by AS^^ < S < AS*^^ is the "intermediate regime" where the system 
alternates between collisional and collisionless reconnection. On the one hand, secondary 
current sheets can reach the di scale and trigger onset of Hall reconnection. On the other 
hand, after onset of Hall reconnection the system cannot settle down to a localized Hall 
reconnection geometry because it is physically unrealizable. This regime has been realized 
in a large scale Hall MHD simulation, where the current sheet becomes extended again after 
onset of Hall reconnection. That leads to formation of new plasmoids and another onset of 



Hall reconnection. 



321 



From what we have learned so far, among these phases, the collisional plasmoid- 
dominated regime gives the slowest reconnection rate, which is approximately O.OIVaB. 
The collisional Sweet-Parker regime yields faster reconnection, because resistivity is high. 
Collisionless Hall reconnection rate typically attains a maximum value ~ O.IVa-B, whereas 
in the intermediate regime the reconnection rate can oscillate between O.OlV^i? and O.IVaB. 



There have been different renditions of phase diagrams proposed in recent literature. 



32 



551-157] All of them share a great similarity, with some minor differences. Because the 
parameter space has not been systematically explored, these phase diagrams are necessarily 
speculative to some extent. Nonetheless, they can serve as a good frame of reference to 



explore large scale magnetic reconnection, either in choosing simulation parameters 



in planning for future experiments. 



321 or 



551] One should bear in mind that the border lines between 



different phases are not ironclad, because there is no clear-cut value of Sc- Furthermore, 
secondary current sheet widths can deviate from the estimated value S ~ LSl^'^/S. Finally, 
both Sc and 6 can be affected by Hall effects. joOl More caveats in applying the phase diagrams 



are discussed in Ref. 



57| . and a compilation of parameters for various astrophysical, space. 



and laboratory plasmas can be found in Ref. 



55|. 
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VI. DISCUSSION AND FUTURE CHALLENGES 



Although significant advances have been made in recent years on large scale, high- 
Lundquist-number reconnection, where plasmoids play important roles, many questions re- 
main open. The results presented here are based on 2D studies. To what extent do these 
2D results carry over to 30 geometry, where oblique tearing modes have been shown to play 
an important role? IgiI. 



62 



turbulent reconnection 



The relationship between plasmoid-dominated reconnection and 



63M66| is another issue that needs further investigation. What are 



their similarities and differences? Will the interaction between overlapping oblique modes 
in 3D lead to self-generated turbulence and further blur the line between the two? 

How do global conditions affect the plasmoid instability and magnetic reconnection in 
general also needs further assessment. One of the most important results from recent 2D 
studies is that reconnection rate cannot be significantly slower than 10~^Va-B, therefore is 
always fast (although this has only been tested up to S* ~ 10^). This conclusion poses new 



57 



67|-l70| 



challenges to theories of solar flares that assume the existence of slow reconnection 
However, these studies were done in simple 2D configurations, and it is not clear whether the 
results can be directly ^applied to solar corona, where the line-tying effect is thought to play 



an important role. 



75 



7lN74| Line-tying is known to have stabilizing effect on tearing instability 



77H79| Both effects could significantly affect 



76| and smoothing effect on current sheets, 
the plasmoid instability. Although line-tying has been employed in the lower boundary of 
some recent 2D simulations of the plasmoid instability in coronal current sheets, 80 



these studies are limited to anti-parallel reconnection, where line-tying stabilization is not 
effective. Line-tying stabilization of the tearing instability is more effective for component 
reconnection, with the guide field line-tied to the boundary. Therefore, to study the effects 
of line-tying on the plasmoid instability, 3D simulations will be needed. 

Even in 2D systems, the current understanding is far from complete. The parameter 
space needs to be systematically explored with higher Lundquist number and larger system 
size than that have been achieved, to test the proposed phase diagrams. The kinetic model 
of plasmoid distribution can be further improved by including other coarse-graining variables 
in addition to ip and f , and the distribution of plasmoids in collisionless regime needs to be 
further studied. 
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